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SIMULATION OF THE GREEN'S FUNCTION 

Kenneth J. Baumelster 


National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

A time-dependent finite difference formulation 
to the Inhomogeneous wave equation Is derived for 
plane wave propagation with harmonic noise sources. 
The difference equation and boundary conditions are 
developed along with the techniques to simulate the 
Dirac delta function associated with a concentrated 
noise source. Example calculations are presented 
for the Green's function and distributed noise 
sources. For the example considered, the desired 
Fourier transformed acousltlc pressures are deter- 
mined from the transient pressures by use of a 
ramping function and an Integration technique, 
both of which eliminates the nonharmonic pressure 
associated with the Initial transient. 

Nomenclature 
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*s 

ax 
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aroa, m? 
amplitude 

range of source, eq. (41) 
function of w 

ramping coefficient, eq. (33) 
dimensionless speed of sound, c*/L*f* 
speed of sound, m/sec 
dimensionless force F*/( p * 0 f*c*) 
dimensionless component of force In x 
direction 
function 
frequency, Hz 

time dependent Green function 
steady-state Green function 
unit step function 
Index of last grid point 
total number of time increments per 
cycle of base frequency 

time Index 

characteristic length, m 
Fourier transform operator 
harmonic number 

unit outward normal at boundary 
time-dependent dimensionless acoustic 
pressure 

spatial ly dependent steady-state 
acoustic pressure 

dimensionless mass source, Q*/f*p* 0 
dimensionless time dependent acoustic 
source, eq. (6) 

spatially dependent acoustic source, 
eq. (9) 

period, sec (1/f*) 
dimensionless time, t*/T* 
time, sec 
time Increment 

maximum stable time Increment, eq. (24) 

acoustic velocity, U*/c* 

space coordinate, x*/L* 

location of source 

space Increment 

(c at/ax)? 

parameter, eq. (A3) 

Dirac Delta function 
exit Impedance, eq. (15) 


e* 0 ambient density, kg/m 3 

o (c at/ax)/ce 

w dlmesloniess frequency (n*>*/f*) 

«*1 2 » 

*n 

dimensionless angular frequency, 
rad /sec 


Superscripts 

* dimensional quality 

k time Index 

Subscripts 

1 space Index 

n harmonic number 

Introduction 


Finite difference and finite element numerical 
solutions to the homogeneous wave equation have 
been successfully used to predict harmonic radia- 
tion from turbofan engine Inlets!*.?) as well as 
Internal propagation In a wide variety of complex 
hard and soft wall ducts, as cited in the summaries 
of references 3 and 4. In these solutions, the 
noise sources were decoupled from the propagation 
problem which reduced the more general Inhomogen- 
eous wave equation' re ** eq. I. 20) t 0 a s t m pi er 
homogeneous form. In the homogeneous form, a known 
or estimated pressure profile is used as a boundary 
condition at some location, such as the fan face 
in a turbofan engine, and the propagation of this 
pressure wave Is predicted from a solution of the 
wave equation or the more general linearized gas 
dynamic equations. 

Numerical solutions to the Inhomogeneous 
equations could be useful in studying the coupled 
Interaction of the noise source and sound field 
for a wide class of practical problems. For ex- 
ample, the noise generated bv propellers' 6 ', fan 
blades'^), and shear flow'®*") can be modeled 
by solutions to the inhomogeneous wave equation 
utilizing monopoles, dipoles, and quadrupoles as 
the noise sources. As with the numerical approach 
to the homogeneous equations, numerical solutions 
of the Inhomogeneous equations would focus on sound 
propagation with complex geometries and flow pat- 
terns''''***) which would be difficult to handle 
with analytical methods. 

Some Inhomogeneous problems of current In- 
terest, such as a turboprop noise field, require 
three dimensional solutions to the sound field. 
Consequently, the transient finite difference 
solution to the wave equation will be employed In 
the present paper. Customarily, In the bulk of 
the past numerical studies' 3 **), the pressure 
and acoustic velocities were assumed to be simple 
harmonic functions of time. In this case, the 
time independent wave equations can be used. The 
matrices associated with the numerical solutions 
to the time independent equation (called steady- 
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state solutions herein) must be solved exactly 
using such methods as Gauss elimination. As a 
result, large arrays of matrix elements must be 
stored, which prohibits realistic three dimensional 
solutions to the wave equation. In contrast, the 
transient analysis begins with a harmonic noise 
source radiating Into an Initially qulesent envi- 
ronment. Next, an explicit Iteration method calcu- 
lates stepwise in real time to obtain the transient 
as well as the steady-state solution of the acou- 
stic field. Advantageously, this approach has 
been chosen because the matrix storage requirements 
are completely eliminated in the transient solution. 

Maestrello, Bayless, and Turkel(^) have 
studied transient finite difference solutions to 
the linearized continuity and momentum equations 
for a pulse noise source which propagates In a jet 
(sheared flow). They found amplification of the 
noise source due to shear. Up to the present time, 
however, transient finite difference solutions 
have not been applied to the harmonic (e“'" *1**) 
form of the Inhomogeneous wave equation. There- 
fore, the primary purpose of this paper Is to 
develop the Inhomogeneous tranlsent finite differ- 
ence equations In one dimension for a harmonic 
sound source with the Somnerfeld radiation boundary 
condition!^) ^ and to verify the theory with 
appropriate examples. 

The transient numerical procedures employed 
herein are intended to simulate the analytical 
Fourier transformed form of the acoustic solutions, 
where all the acoustic parameters vary as e -1 “*l t *. 
The steady-state numerical solutions exactly simu- 
late the analytical results. Unfortunately, this 
is not the case for the transient solution. Basic- 
ally, the analytical and steady-state numerical 
theories assume the harmonic forcing function has 
begun at time equal to minus infinity. On the 
other hand, tn»: transient solution abruptly begins 
the harmonic oscillation at time equal to zero. 
Clearly, 

[ -i«*t*] f -i«.n*l 

e n J*L f [H(t*)e n J (1) 

where Lf Is the Fourier transform operator( ref * 

P* 28) an d h represents the unit step function. 
Symbols are defined In the Nomenclature. The sym- 
bols with an asterisk represent dimensional quan- 
tities. The second purpose of this paper is to 
develop the steps necessary to adapt the transient 
solution to simulate the Fourier transform form of 
the solution. 

Finally, analytical predictions for noise 
propagation from both steady and unsteady noise 
sources, such as a nropeller, rely on the Integra- 
tion of a unit source represented by a free field 
Green's function g and the product of the magni- 
tude of the source F*' re '* P* SuC h that 

P* - J A * F*g*dA* (2) 

Consequently, predictions of the free field Green's 
function as well as the effect of a distributed 
source are developed herein In two examples for 
code verification and assessment of numerical ac- 
curacy and stability. 


Governing Equations and Boundary Conditions 

The governing differential equation and bound- 
ary conditions are Introduced In this section of 
the report. First, the scalar Inhomogeneous wave 
equation Is developed. Next, the Green function 
form of the equation Is developed and will be used 
for code verification. Finally, an expression for 
the Sommerfleld exit boundary condition Is pre- 
sented. 


Inhomogeneous One Dimensional Wave Equation 

In the absence of mean flow, the dimensionless 
form of the plane wave linearized gas equations can 
be written as( r< *- «*• 

Continuity (isentroplc) 




Q 


Momentum 




F 


x 


(3) 


(4) 


The parameters used to reduce the pressures and 
velocities to their dimensionless form are given 
In the Nomenclature. The usual notation for acou- 
stic pressure and velocity are used. Here, Q 
represents a dimensionless mass source and F x 
a dimensionless body force. Differentiating equa- 
tion (3) with respect to time and equation (4) 
with respect to x, and combining yields the di- 
mensionless wave equation 


1 a 2 P 

7 77 


» 2 p 

77 


s (x,t) 


(5) 


where the acoustic source Is 


S(*.t) (6) 

The first term on the right hand side of equation 
(6) represents the monopole contribution while the 
second term represents the dipole contribution. 

The source will be assumed to be harmonic In nature 
(e-in«*lt*) such that equation (5) can be wrlt- 


, 2 P -V 

U- s(x) e ^ H(t) (7) 

»x‘ 


( 8 ) 

first harmonic, n - 1, Is considered 
However, In the numerical data reduction 
a procedure will be presented to handle a 
frequency source problem. Because of the 
Introduction of complex notation for the noise 
source, all the dependent variables are complex. 

Green Function Equation 

For the purpose of checking the numerical 
theory, the simplest form of the source Is desired. 


ten as 

1 a 2 P 

T — 7 ~ 
C ar 

where 

■ n2» 


Only the 
herein, 
section, 
multiple 
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Consequently, s(x) Is assumed to be of the form 
s(x) - s(x s )«(x - x 5 ) (9) 

and 

s(*s) - 1 ( 10 ) 

Thus, equation (7) can be written as 

1 . 2 p - 1 w*t 

— jr — 7 n - g (x - x ) e H(t) ( 11 ) 

c »x' s 

This source Is analogous to the analytical form 
treated by Morse and Ingard' re '* P* *33) f or 
which a closed form analytical solution Is pre- 
sented. 


locltles are taken to be zero. For times greater 
than zero, the application of the noise source 
(eq. (7)) will drive the pressures. 

Difference Equations 

Instead of a continuous solution for pressure 
In space and time, the finite-difference approxi- 
mations will determine the pressure at Isolated 
grid points In space as shown In Fig. 1 and at a 
discrete time steps at. Starting from the known 
Initial conditions at t - 0 and the boundary con- 
ditions, the finite-difference algorithm will 
march-out the solution to later times. 

The difference equations can be developed by 
an Integration of the governing wave equation, 
equation (7), In the form 


If the signal e'M* Is assumed to be ap- 
plied from time minus Infinity, the time dependent 
pressure In equation (11) represents the Green 
function solution for a concentrated source. Thus, 
equation ( 11 ) can be written as 

t ( X _ x ) e " 1 ”! 1 ( 12 ) 

c c »t ax' 5 


t+at/ 2 rx+ax /2 / ^ j2p # 2 p) 


% / 




t-at /2 x-ax /2 

^t+at/ 2 ^.x+ax /2 

t-at /2 x-ax /2 


— 7 7 

ax' 


dxdt 


s(x) e - *"** dxdt 


(17) 


In this steady-state case, all the acoustic para- 
meters vary as e - ’"^, that Is 

6 (x.t) - g(x) e- 1 -^ (13) 

and the differential equation ( 8 ) becomes 


The Integration procedure Is fully documented In 
reference 16 and Is especially useful In developing 
equations at boundaries and for the special case 
of the unit Impulse function. 

Central Cell 



Later, numerical solutions of equation (11) 
will be compared to analytical solutions of equa- 
tion (14). 

S pinner fe Id Radiation Boundary Condition 

The two required boundary conditions upstream 
and downstream of the noise source can be expressed 
in terms of a specific acoustic impedance defined 
as 


Away from the boundaries. In the central cell 
of Fig. 1, the integration of the various terms In 
equation (17) yields 


ax 

cat 


,k*l ?p k + p k-l\ at / p k 
1 ' ' H i *i J ~ ax 1*1 

|r 

- 1 u,t K 

« s(x. )e atax 


"I * p li) 
(18) 


where 1 denotes the space Index, k the time 
Index, and ax and at the space and time mesh 
spacing, respectively. All the spacings are as- 
sumed constant. The time t* is defined as 
t* ■ kat. 



(15) 


Substituting equation (15) Into equation (4) yields 
the following relationship for the exit pressure 
gradient 


if 1 »P 
ax “ " cc e Tf 


(16) 


For the plane wave propagation to be considered 
herein, c e - 1 (dimensionally equal to p<yc) 
which Is equivalent to propagation at an exit with- 
out acoustic reflections, :onmonly called the Som- 
merfeld radiation conditlor. 

Initial Conditions 


In the numerical calculations, for times equal 
to or less than zero, the acoustic field Is assumed 
quiescent, that is, the acoustic pressure and ve- 


Solvlng equation (18) for the new pressure 
Pif 1 yields 

p^ 1 » 2P 1 ; - p 1 ;- 1 - zp 1 ; ♦ 

-lujt _ 2 

♦ Stx^e 1 c'at' (19) 

where a - (cat/ax)?. As will be discussed 
shortly In conjunction with equation (24), the 
value of a must be chosen less than or equal to 
1 for numerical stability. Equation (19) Is an 
algorithm which permits marching out solutions 
from known values of pressure at times associated 
with k and k-1. The procedure Is explldit 
since all the past values of P k and P*“l 
are known as the new values of P* c ' f l are com- 
puted. For the special case at t - 0, the dif- 
ference equation for the start up Is 
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i .0 »p( « /„o „„0 / „0 \ determine the effect of the boundary conditions on 

P 1 * P 1 * 4t lf|o * f ( P 1+l ~ 2P 1 + P 1-l) stability. 


s(x 1 )c 2 at 2 

2 


(20) 


where PQ and »P/»t| 0 are taken as zero, 
and erM* Is taken as 1 and H(t) is 
assumed unity over the entire at 12 time step. 


Exit Boundary Cell 

The solution of equation (17) for the exit 
cells Is fully documented”®) when s(x) Is zero 
(homogeneous solution). In this case, the algo- 
rithm for the exit pressure Is 



2Pj - Pj" 1 + oP 


- 1 -! 


♦ S(xj)c 2 at 2 e 



k 

I 


o) 



( 21 ) 


Numerical Data Reduction 


Recall, at the start of the numerical calcu- 
lation, the acoustic pressures and velocities were 
assumed zero throughout the duct and a source be- 
gins a harmonic oscillation. The pressures will 
vary harmonically at all points In the duct. In 
order to conveniently Interpret the numerical re- 
sults, the time harmonic component of pressure 
should be removed so that only the spatial or 
Fourier transformed component of pressure remains. 

Homogeneous Equation 

Conslde- the solution for plane wave propaga- 
tion In the absence of noise sources (S(x.t) ■ 0) 
and when the pressure at a boundary (x - 0) Is 
given by 

P(0,t) - e-M 1 - < t < - (25) 


where a equals (cat/axe). A similar equation 
exists at the entrance, In this case, I goes to 1 
and 1-1 becomes 2. 

6reen‘s Function Equation 

Assuming s(x) of the form given by equa- 
tions (9) and (10) and noting that H(t) Is 
Identical to 1, the Integration of the source term 
In equation (17) yields 


In this steady-state case, all the acoustic para- 
meters are proportional to e^W* and equa- 
tion (7) reduces to the Helmholtz equation. The 
solution for pressure Is (see Appendix A, Problem 1) 

P(x,t) - e-M(t-*/c) (26) 

However, In the numerical simulation to this exam- 
ple, the boundary condition given In equation (25) 
becomes 


C t-*-at/2/*x+ax/2 -1w,t 

J J s(x)e 1 dxdt 

t-at/2 x-ax/2 

-lo-t* /■x^ax/2 . , t k 

- e 1 at/ a(x-x c )dx - e 1 it 

1-6x12 s (22) 


Therefore, by comparing the results of the Inte- 
gration In equation (22) to the right hand side of 
equation (18), If s(x^) Is chosen such that 


P(0,t) - e-M® H(t) (27) 

In this case, the numerical solution will approxi- 
mate (see Appendix A, Problem 2) 

P(x,t) - e -1-l(t-x/c) H (t-x/c) (28) 

Although only a specific example has beer, con- 
sidered, In homogeneous problems with forcing bound- 
ary conditions, such as equations (25) or (27), the 
pressure P(x,t) will always be proportional to 
e - ’"!® provided sufficient time has elapsed 


s(x^) - 1/iX 


(23) so that the Initial transient has passed (t > x/c). 
Therefore, the numerical data can be reduc 


then equation (19) will simulate the Green's function simple spatial form by dividing by e - '*!® 
solution. Of course, s(x) equals zero for all 1 that is 
values not associated with the position of the con- 
centrated source. 


.If Me ! 0 


St.MHty 


*>»' c 


«9) 


In the explicit time marching approach used 
here, round-off errors cai grow in an unbounded 
fashion and destroy the solution If the time In- 
crement it Is too large or the Iteration schene 
Is Improperly posed. The von Newman methcd( re ** 15, 
P* 104) applied to the homogeneous form of equa- 
tion (7) required that 

itj < ax/C (24) 

This result was first derived by Courant, Fried- 
richs, and Levy In 1927 for a related method and, 
therefore. Is known as the CFL condition. This 
criteria proved successful In the solution of 
equation (7). Note, for a two dimensional problem 
for which nonpropagating modes exist, the exit 
Impedance can Introduce an Instability Into the 
equation(l') . In this case the eigenvalue ap- 
proach( re *- p. 263) should be employed to 


Inhomogeneous Problems 

For many sources, such as a pure monopole, 
numerical data reduction for the Inhomogeneous 
wave equation could employ equation (29) because 
the results will be solely harmonic In time, 
similar In form to equation (25) and (28). How- 
ever, this will not always be the case as will now 
be shown. Consequently the data reduction scheme 
employed must always be scrutinized. 

For the example source considered In this 
paper, the analytical solution to the steady 
Green's function equation (14) Is (see Appendix A, 
Problem 3) 


9(x) 


p(x) 



1c p 1 -l' x - x s' ,c 


(30) 


4 



while the numerical simulation employing equa- 
tion (11) yields (See Appendix A, Problem 4) 


P(*.t) » 

*1 

Dividing by 
yields 


-1*,t lihU-X i/c , 
e -1 e -1 S _ 1C_ t > x/c (31) 

Z “1 

e"M l according to equation (29) 



1c S'*-V /c 


ic p 

‘*=T e 


1*ht 


(32) 


As mentioned In the Introduction, It Is desirable 
that the numerical solutions should approximate 
equation (30). Therefore procedures are developed 
to eliminate the dc component contribution in 
the numerical prediction of p. In this example, 
two new schemes for data reduction will now be 
employed to simulate equation (30). 

Ramping Solution 

The dc component In equation (31) could 
result from the sudden onset of the source at time 
t - 0, In contrast to the Fourier transform solu- 
tion for which the source has been applied from a 
time of minus Infinity. To reduce the dc com- 
ponent, an exponential ramp function (1 - e -bt ) 
has been applied to the source such that 
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ft+i i-.t 

P (*) - J KM) e n dt (35) 


Besides yielding the Individual higher harmonics, 
the orthogonal properties of equation (35) will 
eliminate the steady offset Involved In the numer- 
ical example considered herein. 

In employing equation (35), for convenience, 
the at from equation (28) Is adjusted as follows 
to obtain an even number of at's In one period: 

“ FLOW (inX(i.M s ) riT (36) 

where IFIX converts a floating point number to 
an Interger, and where 

I cyc i e - IFIX (l./atj) ♦ 1 (37) 

Consequently, the following algorithm was used 
to numerically evaluate equation (35) to calculate 
p(x) for n • 1. 



Sample Problems - Comparisons 


S(x,t) - »(x-xj)(l - H(t) (33) 

In this case, the numerical simulation for the 
Green's function solution becomes (see Appendix A, 
Problem 5) 


P(«) 




4 s - e 

Zw, 


i*y x- 


X $ ' Ic 


Two examples are now presented to Illustrate 
the numerical techniques. Both examples have an- 
alytical solutions; therefore, a direct comparison 
between the numerical and analytical techniques 
can be made. The first example presents a simula- 
tion of the Green's function. The second example 
presents a solution for a distributed noise source 

Example 1 - Green's Function 



Equation (34) will approximate the analytical 
Green's function solution provided (b lte\) Is 
chosen to be small compared to unity. 

Integral Solution 

For a large class of Inhomogeneous problems, 
the use of equation (29) presents some additional 
dlflcultles. Consider, for example, a typical 
forcing function for the steady or unsteady (cir- 
cumferential variations of steady /low) loading 
noises on a propeller source' ref * p * 73 °). 

The propeller source can be conveniently modeled 
by stationary sources which turn on and off 
(square wave) with t1me^ re '- *’9* “*19) to 
match the passing of the propeller. For this 
square wave source function In time, the analytical 
expression for pressure will contain all the har-^ 
monies (n«*) of the blade passing frequency e~' n '* * . 
The amplitudes of the various Fourier components 
are determined by an exponential Fourier series fit 
of t^ yi^n l° r »ss»**<l) forcing funct1onl r * f * 

To reduce numerical data which contain a 
spectrum of discrete harmonics resulting from the 
Fourier fit of the square wave, the following 
orthogonality condition is used. 


As our first example of the difference tech- 
nique, consider the numerical solution of the 
Green’s Function equation, (12). In this case, 
the finite difference equations (19) and (20) will 
be used In the solution with the source term equal 
to l./ax (eq. (23)) at x* and zero elsewhere. 

The exact analytical solution Is given by equation 
(32). Taking the position of the source x s 
equal to zero, considering only positive values 
of x In the output, and dividing both sides of 
equation (32) by e“'*lb , equation (32) can be 
rewritten as (recall wj - 2t) 



-Sin 4 1 C0S (“t“) - e1Zl,t 

t > x/c 


(39) 


For a nondlmens lonal speed of sound of unity 
and ax of .00521, the numerical and analytical 
results are shown In Fig. 2 for times of C 75, 
l.G, 1.25, and 1.5. At t - 0.75 in Fig. 2a, the 
Initial transient has progressed three-quarters of 
the way to the termination point of x equal to 
unity. At t greater or equal to 1.0 In Figs. 2b, 
C, and d, the time dependent oscillatory nature of 
the transient form of the solution Is shown. As 
seen In all the Fig. 2 plots, the analytical and 
numerical results are In excellent agreement. 

Recall, the steady-state Fourier transform 
solution to the Green's function equation Is de- 
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sired. In this steady-state solution, the source 
has been assumed turned on from time to minus In- 
finity. The exact analytical solution, eq. (30) 
can be rewritten as 


The results of the Integration technique, 
equation (38) are shown in Fig. 5. As seen In 
Fig. 5, the analytical and numerical results are 
In exact agreement at the end of four cycles 
(t - 4). 


P(x) 



-sin 2*x ♦ 1 cos 2tx 


(40) 


Again, the normalized velocity of sound has been 
taken as unity. 


As discussed in the previous section, the 
numerical aporoxlmatlon to this analytical solution 
can be obtained by either Introducing a ramp func- 
tion or usinq the Integration technique defined by 
equation (35). 


Figure 3 shows the numerical results for the 
exponential ramp with a coefficient of b equal 
to 0.5. In this case about seven time periods 
(t - 7) were required before the transient solu- 
tion associated with the ramp died out and the 
steady-state Fourier transformed solution was 
obtained. For t greater than 7, the numerical 
solution remained unchanged. The analyses, 
equation (34), predicted an error of b/»*i of .08 
which Is the error seen in Fig. 3. 

Finally, the results of the Integration tech- 
nique, equation (38), are shown In Fig. 4. As 
seen In Fig. 4, the analytical and numerical re- 
sults are In exact agreement at the end of two- 
cycles (t - 2). 


In employing the ramping technique for the 
distributed source. It was necessary to reduce b 
to 0.1 to keep an eight percent error, as was used 
In Fig. 3. The reduced b leads to a solution 
time an order of magnitude greater than the Inte- 
gration technique, which limits the usefulness of 
the ramping technique. 

Concluding Remarks 

The transient difference technique was suc- 
cessfully formulated for the one dimensional In- 
homogeneous wave equation with a harmonic source. 
The Green's function and a distributed noise source 
are accurately predicted by the numerical theory. 
The desired Fourier transformed acoustic pressures 
are determined from the transient pressures by use 
of a ramping function and an Integration technique 
both of which eliminate the nonharmonic pressures 
associated with the start-up of the particular 
inhomogeneous source function considered herein. 

The Integration technique is an order of magnitude 
faster than the ramping technique; consequently, 
future .lumerical solution should employ this tech- 
nique. 


Appendix A 

Fourier Integral Solutions 


Example 2 - Distributed Source 

For the case of a distributed sound source, 
the solution for p(x) can be represented by the 
integral ( re *- 14, eq. 4.4.13) 0 f th e product of 
the source magnitude A(x s ) and the Green's 
Function as given by equation (30). In the 
present example, the source strength A Is 
assumed constant between -a and +a and zero 
outside this range. Thus for x greater than a. 


To verify the accuracy and applicability of 
the numerical techniques, closed form analytical 
solutions will now be developed for a variety of 
check problems. The complex exponential form of 
the Fourier integral representation will be used 
to develop the analytical expression. 

The Fourier transform of P(x,t) Is defined 
as (ref. 14, p. 28) 


p(*) - e 
? “1 


1»,x/c 


r' 1 - y 

-a 


a —1 b, x /c 
e 15 dx. 


(41) 


while for x greater than -a and less than a. 


p(x,w) m Yt f e1ut dt (Al) 

— 

The Fourier transform of the wave equation (5) 
becomes 






S(x,t) e^“* dt 


(A?) 


-iw,x/c 
♦ e 


ya e 1 Vs /c 
x 



(42) 


In the specific example now considered, a Is 
set at a value of 0.5, the amplitude A Is taken 
as unity and the dimensionless speed ef sewnd Is 
assumed to be 0.5. In this case, performing the 
integrations In equations (41) and (42) yields the 
following expressions for the normalized pressure 


The Fourier transorms to be considered have 
only simple poles and vanish at t*i ► «*. For 
the problems considered herein, the x and t 
variables can be Incorporatied In a new variable 
$ - t - x/c and the Inverse transform can be 
written as 

i 

P(x,B) - / B(«)e" 1 “ $ d« (A3) 


(4») 2 p ( x ) - -(1 - cos 4» x) x < 0.5 (43) 

p(x) - 0 x > 0.5 (44) 

The parameters have been chosen so that the noise 
sources combine to give complete cancellation In 
the region outside the source. 


I n^ this t|iy be written directly 

S < 0 P ( x , b ) - 2«i I (residues above real 

iu axis) ( A4) 
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6 > 0 P(x,#) - -2*1 l (residues on and below 

real «i axis) (A5) 

Problem 1 - Harmonic Pressure All Times, -•<♦<• 

For s(x) - 0, the solution of equation (A2) 

In a semi-infinite domain (no reflections) yields 


ORIGINAL PAGE 1$ 

OF POOR QUALITY 


p(x,«) - C + e 


1«x/c 


(A6) 


The constant C+ (forward waves) Is determined by 
the transformation of the boundary condition, 
equation (25)C ref - I?. P- 321). 


p(0,n) > a (m - «i) • C+ 


(A7) 


The Inverse P(x,t) of equation (A6) with C+ 
as given by equation (A7) can be determined direc- 
tly from the definition of the Dirac delta func- 
tion^- P- 31* e P- 1-3.24) t0 y1eld 
equation (26) In the body of this report. 

Problem 2 - Harmonic Pressure t > 0 

Again for s(x) ■ 0, the solution of equation 
(A2) in a semi-infinite domain Is given by equa- 
tion (A6). In this case however, the constant 
C+ Is determined by the transformation of the 
boundary condition equation (27) to give 


P(0 


• w) ■ c * 


(A8) 


Problem 4 - Modified Green's Function t > 0 

For a source in the form as given In equation 

iiVMSW S.’V'eTf^r 2 ' 


A.. 2 . 

77 7 0 - TTT.--^T 


(All) 


Using the approach outlined In Problem 3 the 
Fourier transformed pressure becomes: 


p(x,-) 


lull X-X I / C 

i S (A12) 


In taking the inverse of equation (A12), two resi- 
dues are at 0 and and both lie on the 
real axis. Therefore, the inverse can again be 
found directly from equations (A4) and (A5) to 
yield P(x,t) as given by equation (31) In the 
body of this report. 

Problem 5 - Ramping Function t > 0 

For a source in the form as given In equation 
(33) the right hand side of equation (A2) takes on 
the form 

+ . s(x . x $ )[^ r - ~ p- - 77Cr: if* ( b) ] 


Considering equations (A6) and (A8), the only 
residue Is «i, which lies on the real axis. 
Therefore, the Inverse can be found directly from 
equations (A4) and (A5) to yield P(x,t) as given 
by equation (28) In the body of this report. 

Problem 3 - Green's Function - • < t < - 

For a source in the form as given in equation 
(12), the right hand side of equation (A2) takes 
on the formi'' e ^- P- 321 ) 


(A13) 


The solution for the Fourier transformed 
pressure becomes 


p(x,u) 


“ 2(?» ) ^ ti)(« - «i )] 

c [ 1, 1 

2(2*) (in - + 1 b ) J 


-iwi x-x^t / c 


-iui x-x i / c 


(A14) 


— - ^ 4 — j- p ■ -«(x - x $ ) a (w - (d^ ) (A9) 

The homogeneous solution for nonreflected waves in 
thi ^^rticular example can be written as( ref - 34, 

p(x,u) « Ae^f^-xs) /c x > x s 

- Ae-l“(x-x s ) /c x < x$ (A10) 

for a pressure continuous across x equals x s . 
There is, of course, a discontinuity in the slope 
of p(x, w ) at x - x s . The constant A and 
therefore the magnitude of the discontinuity takes 
on the value x > to satisfy the right hand 
side of equation (A§) at x - x s . Following the 
procedure outlined in great depth by Morse and 
Ingard' re ^- 34, p. 133) the constant A is 
found to be A - (ic/u) a (»-■]). Aqain, the 
inverse can be easily determined from the direct 
use of the Dirac Delta function( ref - 34, p. 3l, 
eq. 1.3.24) P(x,t) as given by 

equation (30) in the body of this report. 


In taking the inverse of equation (A14), three 
residues are at 0, «i, and mi-ib and all lie 
on or below the real axis. Therefore, the inverse 
can again be found directly from equations (A4) 
and (A5) to yield 
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For long times, the last two terms in equation 
(A15) can be neglected. In addition. If blm\, 
is assumed to be small, equation (A15) can be ap- 
proximated by equation (34) in the body of this 
report. 
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Figure 1. - Finite difference grid. 
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Figure 2. ■ Nature of Green’s function and numerical solutions when source is turned on at time equal to zero. 
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Figure 3 - Green function solution to the inhomogeneous Figure 4. - Green function solution to the inhomogeneous 

wave equation using a ramped source (b-0.5 and t-7.0). wave equation using the integration technique (t-2). 
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